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The Helium Abundance in the Ejecta of U Scorpii 
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ABSTRACT 

U Scorpii is a recurrent nova which has been observed in outburst on 10 occasions, 
most recently in 2010. We present near-infrared and optical spectroscopy of the 2010 
outburst of U Sco. The reddening of U Sco is found to be E{B ~ V) ^ 0.14 ± 0.12, 
consistent with previous determinations, from simultaneous optical and near-IR obser- 
vations. The spectra show the evolution of the line widths and profiles to be consistent 
with previous outbursts. Velocities are found to be up to 14000 kms~^ in broad compo- 
nents and up to 1800 kms~^ in narrow line components, which become visible around 
day 8 due to changes in the optical depth. From the spectra we derive a helium abun- 
dance of 7V(He)/A^(H)= 0.073 ± 0.031 from the most reliable lines available; this is 
lower than most other estimates and indicates that the secondary is not helium-rich, 
as previous studies have suggested. 
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1 INTRODUCTION 

U Scorpii is a recurrent nova (RN) which has under- 
gone recorded outbursts in 1 863, 1906, 191 7, 1936, 1945, 
1969, 1979, 1987, and 1999 l|Schaefeij [2OI0I ) before a fur- 
ther outburst peaked on 2010 January 28.19 ± 0.17 UT 
JSchaefer et al.l I2OI0I} . The mean recurrence time is 10.3 
years I Schaefen 2010l ). U Sco i s an eclipsing binary with 
an orbital inclination of ~ 82° (|Thoroughgood et al.|[200ll ) 
and is semi-d etached with an orbital period of ~1.23 days 
ISchaefer fc R ingwald 1995). U Sco consists of a white dwarf 
(WD) primary and a probable s ubgiant secondary with 
a spe ctral typ e in the ran ge K2 (|Anupama fc DewanganI 
l2000l 'l to GO l|Hanesl 119851 ). with the white dwarf hav- 
ing a mass close to the Ch andrasekhar limit (~ I.37M0) 
l|Thoroughgood et al.|[200ll ). The system is at a distance of 
12 ± 2 kpc and is far out of the galactic plane at a height of 
~ 4. 5 kpc ('SchaefeJ l2oTol ). 
IStar rficld (20Qi) 

interprets the outbursts as being due 
to a thermonuclear runaway (TNR) on the surface of the 
white dwarf. The TNR occurs when the temperature and 
density at the base of the layer accreted from the secondary 
reach critical values of ~ 10* K and ~ 10^® Nm~'^ respec- 
tively (Starrficld 2008 ). The recurrence time scale is consis- 



Email: mpmaxwell@uclan.ac.uk 



tent with the nova outburst models of lYaron et all (|2005h . 
The energy released from the hydrogen burning is suffi- 
cient to allow heavier elements to undergo nuclear fusion. 
This continues until the energy generation is limited by the 
long, temperature independent half-lives of some isotopes 
involved in the CNO cycle such as and ^'^N. 

The 2010 outburst was anticipated bv lSchaefeil (j2004l ) 
who planned a multiwavelength observing programme ahead 
of time; this resulted in the 2010 outburst of U Sco having 
the best temporal coverage of any nova event so far. This 
led to the discovery of new pheno mena such as aperiodi c 
dips in the light curve and flares (|Pagnotta et al.l l|201ll ) ; 
Schaefer et al. in prep). The a pparent onset of op tical flick- 
ering on day 8 of the outburst jWorters et al.|[2010h has been 
identified as an example of such a flare. The spectral evo- 
lution between outbursts is consistent, as is the photomet- 
ric evolution. Bancrjcc ct al. (2010') present near-IR spec- 
tra of the outburst which show broad Hi, He I, Hell, and 
O I emission lines. The H I lines give an upper limit on th e 
ejected mass of 9.71 ±9.29 X 1O"^M0 (|Baneriee et al.ll2010l ). 
The heli um abu ndance, A'^(He)/A'^(H), of U Sco is highly un- 
certain jDi^_et^al. 20103 with estimat es ranging from 0.16 
(|liiimall2002h to 4.5 (|Eva,ns et al.ll200ll ). An accurate deter- 
mination of the helium abundance in U Sco is necessary as 
some studies have suggested that, unlike classical novae, the 
WD is accreting helium-rich material from the secondary 
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star. A s a result, the physics of the TNR has been modi- 
fied bv lStarrfield et al.l l| 19881 ) for the case of U Sco. To ex- 
plain the pres ence in the system of a hehum-rich donor star, 
iHachisu et al.i (,1999 ) proposed an evolutionary sequence in 
which a helium-rich primary filled its Roche lobe and trans- 
fered material to the secondary in the post common envelope 
phase. However, iTruran et al.l l|l988l ) showed that a bright 
nova outburst is only possible on a 1.38 solar mass WD if 
H/He < 1. This agrees with evolutionary calculations by 
ISarna et al.l (|2006l ). who concluded that any helium enrich- 
ment is from the WD. 

Here we present near-infrared (IR) and optical spec- 
troscopy of the latest outburst obtained at the New Technol- 
ogy T elescope (NTT), Liverpool Telescope (LT. ISteele et all 
(|2004l )). Cerro Tololo Inter- American Observatory (CTIO), 
and South African Astronomical Observatory (SAAO) from 
which we derive the reddening, helium abundance, and 
broad and narrow component line widths. 



2 OBSERVATIONS 

We take the observed maximum at 2010 January 28.19 to 
be day of the outburst and all epochs are relative to this 
time. Near-IR spectroscopy of U Sco was obtained on days 
5.41 and 9.43 at the NTT using the Son Of ISAAC (SOFI) 
spectrogra ph. The spectral reso lution is ~ 1000. For more 
details see iBaneriee et aP (|2010h . An optical spectrum was 
also obtained on day 9.43 using the Cassegrain spectrograph 
on the 1.5 metre SMARTS Telescope at CTIO at a resolution 
of ~1500. Optical spectra were obtained from 1.93 to 11.93 
days after outburst using the Cassegrain spectrograph on the 
SAAO 1.9 metre telescope with spectral resolution ^1000. 
Optical spectra were also obtained fr om 6.81 to 12.81 days 
after outburst using the FRODOSpec (|Morales-Rueda et al.l 
l2004h spectrograph at LT at a resolution of 5400. 



3 RESULTS 

Figures 1 and 2 show the dereddened optical spectra of U 
Sco covering most of the first 13 days of the outburst. The 
strong emission lines are due to Hi, He I, Nlll, and, at later 
times. Hen. The Hel lines fade as the Hen lines develop and 
are no longer detectable by day 11.81. Nlll and H7 also fade 
and are very weak or undetectable by this time. H5 and He 11 
4686A are both blended with Nlll. Figures 3 and 4 show 
near-IR spectra of U Sco; they show the Paschen series of 
hydrogen in emission along with O I, He I, and He 11 emission 
lines. Paschen 7 and Paschen 5 are blended with Hel and 
Hell respectively. Figure 3 also shows a spectrum from the 
outburst in 1999 taken at a similar time after maximum to 
the first of our NTT spectra; the spectra are very similar 
with the same emission features present and similar relative 
line strengths. 



Table 1. 


Observing Log 




Day 


Wavelength range (A) 


Facility 


1.93 


35UU-725U 


SAAO 


4.93 


35UU-725U 


bAAO 


5.41 


995U-24UUU 


NTT 


5.93 


35UU-725U 


SAAO 


6.81 


3900-5200, 5700-7900 


LT 


7.81 


3900-5200, 5700-7900 


LT 


7.93 


3500-7250 


SAAO 


8.81 


3900-5200, 5700-7900 


LT 


8.93 


3500-7250 


SAAO 


9.43 


4000-4800, 9950-24000 


CTIO, NTT 


9.93 


3500-7250 


SAAO 


10.93 


3500-7250 


SAAO 


11.81 


3900-5200, 5700-7900 


LT 


11.93 


3500-7250 


SAAO 


12.81 


3900-5200, 5700-7900 


LT 



3.1 Line Widths and Flux Ratios 

Here we measure the FWZI of each line present in our spec- 
tra. Some lines also show narrow components; in these cases 
we also measure the FWHM of that component. The contri- 
bution of instrumental resolution was not found to be signif- 
icant. Velocities and the associated errors were measured by 
fitting Gaussian profiles and visual inspection of where the 
emission meets the continuum. The Hq veloc ity is consis - 
tent with previous outbursts at similar times l|liiima|[2002l ) 
at 8000 — 9000 kms~^. The width of each line can be seen 
in Tables 2-4, and can be seen to be changing with time, 
however the profile of the lines is changing rapidly as can be 
seen in Figures 1 and 2. Tables 3 and 4 show the evolution 
of the narrow components of Ha and H/3. Such high velocity 
outflows (~ 10,000kms"^) were also s een in early (<5days ) 
spectra of the 2010 outburst of U Sco (|Baneriee et al.ir2010l ). 

Line fluxes were measured by fitting Gaussian profiles 
to the lines with a least squares fit and can be seen in Table 
5. Multiple profiles were used in some cases; an example 
fit can be seen in Figure 5. These fluxe s were then ratioed 
and co mpared to theoretical values from lHummer fc StorevI 
From these ratios the abundance of species A to 
species B can be found using the equation 

-Fa _ Aa Na_ Qs 

Fb ' Xb Nb aA' ^ ' 

where F is the measured dereddened flux, A'^ is the abun- 
danc e by number, and a is t he recombination coefficient 
from iHummer fc StorevI l| 19871 ) for species A and B respec- 
tively. We assume that the sources of Hi, Hel, and Hell 
emission are co-extensive. 

This process was repeated for each available Hel (for 
He""") and Hell (for He"""^) to hydrogen line ratio. Since our 
spectra show both ionised states of helium in emission the 
total abundance can be directly measured instead of using 
the Saha equation to predict the abundance of one species 
from the other. 
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Figure 1. Spectra of U Sco taken from day 1.93 to 12.81 at LT (days ending .81), CTIO (day 9.43), and SAAO (days ending .93). The 
spectra are offset as indicated and the flux is in arbitrary units. LT spectra are smootlied witli a Gaussian profile. 



Table 2. Broad Component Line Widths (FWZI, kms-^) from LT (days ending .81) and SAAO (days ending .93). 



H7 N III 4616A He 11 4686A H/3 Ha He i 6678A He i 7065A 



1.93 






8000±2000 


lOOOOitlOOO 






4.93 






9000±2000 


9000±1000 






5.93 


12000±2000 




12000±2000 


9000±1000 






6.81 


9000±3000 


8000±2000 6000±2000 


9000±2000 


8000±1000 


2000±500 


13000±3000 


7.81 




8000±1000 6000±2000 


9000±2000 


8000±1000 


1000±500 


10000±2000 


7.93 


10000±2000 


6000±3000 


9000±2000 


9000±1000 






8.81 


8000±2000 


12000±2000 4000±2000 


8000±2000 


8000±500 


1500±250 


12000±2000 


8.93 


9000±2000 


6000±2000 


9000±1000 


9000±1000 






9.93 


4000±1000 


4000±1000 


3000±1000 


9000±2000 






10.93 


4000±1000 


2000±500 


3000±1000 


9000±1000 






11.81 




1500±500 


2000±1000 


8000±2000 






11.93 


dOOOilOOO 


2000±.500 


2000±.500 


9000±1000 






12.81 




1200±250 


ir)()0±r,()o 


SOOOilOOO 






Table 3. Narrow Component Line Widths (FWHM, kms"!) from LT data. 


Day 


6.81 


7.81 


8.81 


11.81 


12.81 





H/3 1001.95 ± 72.36 836.21 ± 131.3 953.53 ± 41.62 882.85 ± 61.97 1234.63 ± 85.49 
Ha 1071.35 it 50.0 1146.4 ± 46.77 849.94 ± 23.94 753.96 ± 57.1 1380.0 ± 41.95 
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Figure 2. Spectra of U Sco taken from day 1.93 to 12.81 at LT (days ending .81) and SAAO (days ending .93). The spectra are offset 
as indicated and the flux is in arbitrary units. LT spectra are smoothed with a Gaussian profile. 



Table 4. Narrow Component Line Widths (FWHM, kms^^) from SAAO data. 



Day 


1.93 4.93 5.93 7.93 


8.93 


9.93 


10.93 


11.93 


H/3 
Ha 


1782.15±106.9 
1413.54±68.99 


1318.17±108.8 
876.78±40.20 


975.48±71.76 
778.07±32.47 


764.72±53.5 
841.79±30.53 


1192.30±80.41 
882.26±36.28 



3.2 Reddening 

Several studies (iBarlow et al.lll98ll : lAmores fc Lepindlioosl : 



Burstein fc Heile j llQsj T as well as the dust maps of 
Schleeel et al.l ( 19981 ) have shown the reddening of U Sco 
to be in the range E{B -V) = 0.09 - 0.36. The ratios of line 
fluxes were es timated and compa r ed to the theoretical values 
derived from iHummer fc StorevI (|l987l ). We use the optical 
spectra taken on day 8.81 at LT and day 9.43 at CTIO (Figs 
1 and 2), and the IR spectrum taken on day 9.43 at NTT 
(Fig 4) to estimate the reddening, as by this time the flux 
ratios are converging on case B v alues. We use the extinction 
law and assumption of 7? = 3.1 o f lHowarthI l)l983l ). Using the 
ratios H/3/Pae, H/S/Pa/?, H7/Pae, and H7/Pa;g the redden- 
ing was found to be in the range E{B — V) = 0.0 — 0.29 



with a mean of E[B — V) = 0.14 ± 0.12, consistent with the 
previous studies. Although U Sco is at a distance of 12kpc, 
low reddening is consistent with both the line of sight leav- 
ing the plane of the galaxy and the systern being at a height 
oi z = 4.5 kpc above the galactic plane (Sch acfcr 2010). For 
the spectra used in this work we adopt E{B — V) = 0.2 
as a good compromise between our own value and previous 
estimates. 



3.3 Helium Abundance 

The helium abundance of U Sco was calculated using Hi, 
He I, and Hen recombination lines. We use the line fluxes 
on days 8.81 (LT), 9.93-11.93 (SAAO), and 9.43 (NTT) to 
estimate the helium abundance, as by this time the line ra- 
tios are converging on case B values. The large Ha/H/3 ratio 
throughout the time coverage we have available shows that 
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Figure 3. Top: IR spectrum of U Sco taken at the NTT on day 5.41 of 201 outburst. Bottom: IR spectrum of U Sco taken at the NTT 
5.33 days after the maximum of the 1999 outburst from lEvans et aL 1 1 I2OOII) . offset as indicated on right hand side. 

Table 5. Dereddened line fluxes. Fluxes and errors are relative to H/3 on that date for optical data and Pa/3 for near-IR 
data. 





H/3 


H7 


Pa/3 


Pae 


He I 6678A 


He 11 4686A 


He II 5411 A 


He II 1.163Atm 


8.81 


litO.Ol 


1.02±0.03 






0.007±0.001 








9.43 






1±0.04 


0.45±0.02 








0.18±0.01 


9.93 


1±0.04 


0.29±0.02 








0.34±0.02 


0.074±0.010 




10.93 


1±0.02 


0.12±0.01 








0.19±0.01 


0.022±0.002 




11.93 


1±0.03 


0.25±0.04 








0.28±0.03 


0.036±0.006 





Ha should not be used in our abundance analysis due to 
optical depth effects. 

Abundance analyses of helium are complicated by 
the metastability of the lowest triplet level of Hel, 
which can cause some lines to become opti- 
cally thick. Furthermore, collisional excitations from this 
level enhance triplet lines. The collisional contributions 
are 56% for 5876A, 78% f or 7065A, and 72% for 
1.08 3/jm using relations from [Kingdon fc Ferlandl (Il995l ) 
and IPeimbert fc Torres-Peimbert (1987) at a temperature 
of 2 X 10* K at the high density limit. The contribution is 
lower in the singlet lines, although collisions from the 2^S 
level must be taken into account. The ratio of 7065/4471 



is an optical depth indicator (jOsterbrock fc Ferlandl I2OO6I ) 
and this ratio (1.31) is inconsistent with the case B value 
(~ 0.5), further evidence that the 7065A line is unsuitable 
for this analysis. In our abundance analysis we consider the 
singlet He I line at 6678A for wh ich the collisional contribu - 
tion is 27% using equation 12 of lKingdon fc Ferlandl (|l995l ) 
at a temperature of 2 x 10* K. At temperatures of 10"* and 
3 X 10* the collisional contributions for this line are 7% and 
36% respectively. We also use HeH lines at 4687A (LT and 
SAAO), 54IIA (SAAO), and 1.163/xm (NTT). 



Following lEvans et all l|200lD we estimate the electron 
density using Paschen and Brackett series H I lines from our 
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Figure 4. IR spectrum of U Sco taken at the NTT on day 9.43 of 2010 outburst. 



near- 



-IR spectra (Figs 4 and 5) using the equation: 



F{H) 



r2 y^3 
X 



3D2 



(2) 



where F is the flux of a hydrog e n hne , e is the case B emis- 
sivity from iHummer^ fc StorevI (|l987l ). a is the number of 
electrons per hydrogen, is the volume filling factor, A'^e is 
the electron density, V is the velocity, D is the distance, and 
t is th e time since maximum in days. We follow [Evans et al] 
l|200ll ) by assuming q = 1, </> = 0.01 and D = 1 2kpc . 
Our results are consistent with those of lEvans et al.l (|200ll ') 
at similar times, therefore we adopt their values of elec- 
tron density of 10^" cm~^ for data taken before day 9.43 
and 10^ cm~^ for data taken at later times. The recombi- 
nation coefficients are relatively insensitive to temperature; 
we assume Te = 20, 000 K. We derive 7V(He+)/iV(H+) = 
0.012 ± 0.015 and N{Re++) / N[Yi+) = 0.061 ± 0.010 from 
the data obtained at LT, NTT, and CTIO using the He I line 
at 6678A and the Hen line at 1.163^im. From the SAAO data 
we derive Af(He++)/iV(H+) = 0.047 ±0.011 using the Hen 
line at 4686A and iV(He++)/Af(H+) = 0.076 ± 0.023 using 
the line at 541lA. These results are summarised in Table 6. 



Table 6. Helium abundances 



Helium Line 


Derived Abundance 


He I 6678A 


0.012 ± 0.015 


Hen 4686A 


0.047 ±0.011 


Hen 541lA 


0.076 ± 0.023 


Hen 1.163Atm 


0.061 ± 0.010 



4 DISCUSSION 



The helium abundance in U Sco has been the subject 
of many studies which h ave resulted in a wide range of 
values from 0.16 to 4.5 (lliiimal [20o3: iBarlow et al.l Il98ll : 



lAnupama fc DewanganI l2000l : lEvans et al.l l200ll ). We find 
the abundance to be Af(He) / A^(H) = 0.073 ± 0.03 1 lowe r 
than that found by lliiimal (|2002l ) iBarlow et all JlQSj^ , 
lAnupama fc DewanganI l|2000l 'l. and lEvans et all (|20oii rA 
high helium abundance would suggest that the secondary 
in U Sco is a highly evolved star, however the value de- 
rive d here is cl ose to the solar value N(He)/N(H) = 0.085 
(Asplund et al.i 2009). Our estimate suggests that assump- 
tions about the helium-rich nature of the secondary are un- 
founded and that the secondary did not accrete helium-rich 
material significantly in the post common envelope phase 
iHachisu et al.l ri999 ^V In the evolutionary calculations of 
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Figure 5. Fit of H/3 using a pair of Gaussian profiles. 

ISarna et all (l2006l l. U Sco is best represented by their se- 
quence B, with a hydrogen-rich secondary. For our hehum 
abundance analysis we have considered the data (>day 8) 
and the hnes for which we can be confident that case B re- 
combination apphes. Other estimates have been larger than 
our value, especially those based only on the Hell lines. 
These large abundances result from the use of the Saha 
equation, which is highly sensitive to the assumed electron 
temperature. However, we use both He I and Hell lines in 
our analysis, avoiding the use of the equation to derive the 
total helium abundance. The temperature dependencies in 
this analysis are in the recombination coefficients, which are 
relatively insensitive to the value chosen, and in the coUi- 
sional correction for He I lines. Using a range of lines from 
days 8.81 and 9.43 and a range of recombination coefficients 
we determine several He"'"^/He^ ratios, each of which sug- 
gest electron temperatures of 2.3 — 2.4 x 10^ K, therefore we 
use recombination coefficients for the closest temperature 
from lHummer fc StorevI l|l98'<t ) of Te = 2 x 10'' K. 

The hydrogen line profiles have broad and narrow com- 
ponents. The narrow components become prominent around 
day 8 and arise from slow moving material at speeds of 
~ 1000 kms~'. The broad components are composed of ma- 
terial travelling at velocities of up to ~ 10, OOOkms"^. The 
rise in the strength of the narrow components could be ex- 
plained by optical depth effects. The ejecta begins to become 
optically thin around day 8; before this time slower mov- 
ing material would be concealed by the shell of optically 
thick fast moving material. As the faster moving material 
expands and becomes optically thin the flux from the slow 
moving material becomes visible. Hen lines also begin to 
dominate the helium emission around day 8, however this 
is most likely due to photo-ionisation of helium rather than 
an optical depth effect. The measured line widths are gen- 
erally higher for the SAAO data than for near-simultaneous 
LT data; possible reasons for this are the difference in res- 
olution and signal to noise issues. There also appears to be 
clear structure to the H/3 line in the LT spectra which is not 
present in the SAAO spectra. A possible source of the nar- 
row components is an inclination effect; since U Sco is a high 
inclination system most of the velocity is in the plane of the 
sky, therefore we may be observing a small radial compo- 



ne nt. This is consist e nt wi th the asymmetric ejecta models 
of iDrake fc Orlando! (|2010l ). Data over a longer time range 
and with better orbital phase coverage is required to explore 
this further. 



5 CONCLUSION 

We present optical and near-infrared observations of the 
2010 outburst of U Scorpii. We find the helium abundance of 
U Scorpii to be N{Ue)/N{U)= 0.073 ±0.031. This estimate 
is lower than most previous studies and does not support 
their conclusions which suggest that the secondary in this 
system is helium-rich. The velocities seen in the 2010 out- 
burst are consistent with those seen in previous outbursts, 
with some hydrogen and He ll lines being seen to have both 
narrow and broad components in our later spectra. Fur- 
ther observations are planned to investigate the nature and 
metallicity of the secondary in this system. 
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